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Abstract 



The (non-strange) vector and axial-vector spectral functions extracted from r- 
decay by the ALEPH collaboration are confronted with QCD in the framework 
of a Finite Energy QCD sum rule (FESR) involving a polynomial kernel tuned 
to suppress the region beyond the kinematical end point where there is no longer 
data. This effectively allows for a QCD FESR analysis to be performed beyond 
the region of the existing data. Results show excellent agreement between data 
and perturbative QCD in the remarkably wide energy range s = 3 — 10 GeV 2 , 
leaving room for a dimension d =4 vacuum condensate consistent with values 
in the literature. A hypothetical dimension d=2 term in the Operator Product 
Expansion is found to be extremely small, consistent with zero. Fixed Order and 
Contour Improved perturbation theory are used, with both leading to similar 
results within errors. Full consistency is found between vector and axial- vector 
channel results. 
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1 Introduction 



Twenty five years ago [JJ it was pointed out that the hadronic decay of the r- 
lepton constitutes an ideal laboratory for studying the hadronic weak currents at 
low and intermediate energies. It was argued that the inclusive character of this 
decay would make it possible to carry out reliable theoretical calculations in the 
framework of perturbative QCD (PQCD). Vector and axial- vector final states 
would be experimentally identified, thus opening new windows for the study of 
chiral-symmetry. Also, strange and non-strange channels would be separated, 
allowing for a determination of the strange quark mass. Last, but not not least, 
the strong coupling constant would be determined from the total decay rate 
with impressive accuracy. The first experimental results from the ARGUS col- 
laboration at DESY [5] were subsequently used [3] to extract the values of the 
QCD vacuum condensates entering the Operator Product Expansion (OPE) and 
QCD sum rules [1J. When a high statistics r-decay experiment was to become 
a reality at the ALEPH 5j and OPAL [B] detectors at LEP, more detailed stud- 
ies of r-decay were performed [7]. Following these pioneering analyses, literally 
hundreds of papers discussing various theoretical aspects of r-decay have been 
published [8]. One of the key issues is that in the limit of vanishing up- and 
down-quark masses the vector and axial-vector spectral functions should be- 
come asymptotically equal. This is not quite yet the case for the experimental 
ALEPH spectral functions, even at the highest available energies. Therefore, 
PQCD cannot directly be used for accurate calculations in the vicinity of the 
positive real s-axis (s being the energy squared). There is, however, sufficient 
evidence that PQCD and the OPE can be used with impunity for large \s\ in 
the space-like region, as well as in the complex s-plane away from the positive 
real axis. This can be achieved e.g. by using QCD Finite Energy Sum Rules 
(FESR) with so-called pinched kernels which vanish at the end of the spectrum 
[9 - 10J . Some of the highlights of this procedure are (i) a precise determination 
of the strong coupling constant [SJ, [HI-IE], (ii) the precocious saturation of 
the Weinberg sum rules [10], [13], and (iii) the determination of some of the 
vacuum condensates of the OPE [T3J-JT5]. In connection with the latter, a care- 
ful QCD-FESR analysis [15 of the vector and, separately, of the axial-vector 
(non-strange) channel shows that the values of the gluon condensate (a s G 2 ) are 
equal within errors, as expected. This follows, though, only if pinched FESR 
are employed, as otherwise results from the two channels have opposite signs as 
observed in [5]. In addition, as it has been reiterated several times [TO], [T5] , 
[T7] the actual values of the condensates are strongly correlated with the value 
of a s (m 2 ), or alternatively with the value of the QCD scale A. This analysis 
[15] has also found no clear evidence for a dimension d — 2 contribution in the 
OPE, thus confirming results from earlier analyses using the ARGUS data [18], 
as well as the ALEPH data [15], [17], and leading to values consistent with zero. 
The conclusions that the gluon condensate from the vector channel is compat- 
ible with that from the axial-vector channel, and that a dimension-two term is 
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consistent with zero, are less apparent from analyses which do not address the 
convergence problems near the positive real s-axis [5].[19j. 

It is an unfortunate fact that the r lepton is not massive enough to allow for a 
direct confrontation between data and QCD at higher energies. In this paper 
we attempt to overcome this restriction by introducing a new FESR involving 
an integration kernel designed to suppress the hadronic contribution beyond the 
kincmatical end-point of r-decay, si = M%, where there is no longer experimen- 
tal data. We concentrate first on the axial-vector channel, as it involves the 
well known pion pole. Hence, the weighted integral of the hadronic data up to 
s = si, plus that of PQCD extended beyond this point to s = sq > s\, can 
be confronted with the pion decay constant. Results indicate a remarkable and 
stable agreement between QCD and the ALEPH data over the unusually wide 
range so — 3 — 10 GeV 2 , with room for a dimension-four condensate in line 
with independent determinations. A dimension-two term in the OPE, if at all 
present, turns out to be extremely small. The techniques of Fixed Order PQCD 
(FOPT) as well as Contour Improved PQCD (CIPT) [20] have been used, with 
both leading to very similar results within errors. Finally, we study the vector 
channel along the same lines, except that now there is no scalar meson pole 
to set a scale. In spite of this, we find perfect agreement with the conclusions 
drawn from the axial-vector channel. In particular, the value and sign of the 
gluon condensate from this channel is totally consistent with that from the axial- 
vector channel. 



It should be made clear that the purpose of this analysis is definitely not to 
determine the pion decay constant nor the vacuum condensates. The former is 
well known from experiment, and thus supplies a reference scale, and the latter 
have already been determined with different QCD FESR [T3]-[TB]. Instead, our 
purpose is to show how to extend a QCD FESR analysis beyond the kinematical 
region of the existing data. This technique can then be used to confront the 
data with QCD. The conclusion from this analysis is that the ALEPH T-decay 
data in both the vector and the axial-vector (non-strange) channels is perfectly 
consistent with QCD. 



We begin by defining the axial-vector current correlator relevant to r-decay, 



2 The Finite Energy Sum Rule 
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where A ft (x) = d{x)"f tl ^u{x). The spectral function used by the ALEPH col- 
laboration, a^Q (s), is related to the imaginary part of the correlator by 



Imn^( S ) = i-a 1)0 ( S ) 



(2) 



In this paper we will be concerned first with the axial- vector correlator lT^ L+0 ' ) (s) 
which we will simply call II(s) for short. The OPE of n(s) beyond perturbation 
theory can be written as 



4^ 2 n(Q 2 ) = J2 
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Q 
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(3) 



where Q 2 = —q 2 , and q 2 is large and space-like. The term with N = above 
corresponds to PQCD, and the rest of the series parametrize non-perturbative 
physics. In the standard picture the term with N=l is absent from the OPE as 
there is no gauge invariant operator of dimension d = 2, other than negligible 
quark mass terms, that could be constructed from the quark and gluon fields. To 
facilitate comparison with alternative conventions in the literature, absorbing 
the Wilson coefficients (including radiative corrections) into the operators in 
Eq.(3) leads to 
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The spectral function in PQCD to five-loop order, renormalized at a scale fj,, is 
given by 
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is (21] 



0.0862 n^, and the constant £3 



n F ), ft = -|(1Q2- 
0.1153 n F , F± = 

49.076 has 



3 



7T 7T \ TT J \ Pi 



hxL 



+ 



'ai 1} (soA 2 



(In X - InL - 1) + 



ai 1; (s ) 



4 r 



^(ln J L--ln^-21nL+-) 



„/3 2 /3 3l r , &3 
+ ~0f Jh 



(6) 



where 
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with L = ln(so/A 2 ) defines the standard MS scale A to be used here, and 
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with ( 3 = 1.202. Invoking Cauchy's theorem in the complex energy-squared 
s-plane leads to the standard FESR 
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where the dimensionless PQCD moments M 2 jv+2(so) are given by 
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Figure 1: Results for /„■ from the standard FESR in the axial- vector channel, 
Eq. (9) with N=0, with no dimension d = 2 term, and using CIPT, with 
A = 365 MeV (a s (M^) = 0.335). The straight line is the experimental value of 



In Fig.l we show the result for f„ from the standard FESR, Eq.(9) with N=0, us- 
ing the ALEPH data and CIPT with A = 365 MeV [a s (M?) = 0.335), compared 
with the experimental value f n — 92.4 ± 0.1 MeV. Except possibly near the 
kinematical end-point, the agreement is not particularly encouraging. In fact, 
it has been known for quite some time that standard FESR in the axial-vector 
channel do a poor job in reproducing the experimental value of / T . Indeed, 
both the first and the second Weinberg sum rule arc not well saturated by the 
ALEPH data, unless one introduces polynomial integration {pinched) kernels in 
the FESR [H]-[TU]. This may well be an indication that quark-hadron duality 
does not hold on the real s-axis even at higher energies, except probably very 
near the kinematical end-point s\ — M%. Given the absence of experimental 
data beyond this point, we propose the introduction in the FESR of an integra- 
tion kernel, a polynomial P(s), designed to eliminate the (unknown) hadronic 
contribution to the integral between si and some sq > s±. The degree of this 
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polynomial P(s), while in principle arbitrary, should not be too high as each 
additional power of s brings into the FESR an extra higher dimensional con- 
densate with a higher uncertainty. We have found that the optimal degree is in 
fact the simplest, i.e. a linear function 



P(s) 
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so that 
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In this case the complete FESR becomes a linear combination of a dimension- 
two and a dimension-four FESR, which from Eqs.(9) and (11) it is given by 
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where the pion pole has been separated from the data, and the chiral limit is 
understood. The term A(so) is the error being made by assuming that the data 
is constant in the interval Si — s - While the data is unknown in this region, 
if one were to assume the onset of PQCD beyond s = si, then the data would 
follow the logarithmic fall-off of PQCD. In the next section we will address this 
issue and quantify this error, which turns out to be very small. Notice that due 
to the specific form of the integration kernel, P(s), both condensates enter with 
the same overall sign in the FESR. The cutoff point s\ in Eq.(ll) will be chosen 
as close as possible to the kinematical end-point of the spectrum, in order to 
include as much data as possible. The sensitivity of the results to this choice 
will be carefully quantified in the next section. Another point to be considered 
is the fact that both C2{0 2 ) and 64(04} enter in the above FESR. Hence, a 
study of their correlation will be done in the next section. The FESR in the 
vector channel follows from Eq.(13) by simply setting its left hand side equal to 
zero. 
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3 Results 



We begin by evaluating the right hand side of the FESR, Eq.(13), in FOPT 
and comparing the results with the experimental pion decay constant, f n = 
92.4 ±0.1 McV. We have chosen here, and in the sequel, a cutoff point s% — 
2.7125 GeV 2 , and defer to the end of this section a discussion about sensitivity 
to this choice. We recall that in FOPT the strong coupling a s (so) is frozen in 
Cauchy's contour integral, and the Renormalization Group (RG) is implemented 
after integration. In CIPT, instead, the coupling is running and the RG is used 
before integrating. Details of these procedures, as well as explicit formulas 
for the moments Mn(sq) may be found e.g. in 15J. In order to gauge the 
impact of the vacuum condensates, we have first set both condensates in Eq. 
(13) equal to zero and determined / ff using the value of the strong coupling 
obtained in FOPT .lj, i.e. a s (M 2 ) = 0.322 ± 0.004 ± 0.02 (A = 340 - 347 
MeV). Results are shown in Fig. 2. Next, we set 02(02) equal to zero, and 
6*4(04} = 0.10 ± 0.05 GeV 4 , according to earlier independent determinations 
[T5] . [I7]-[T5]. In this case, the results are shown in Fig. 3, which compared with 
Fig. 2 clearly indicates that the contribution of the dimension-four condensate is 
rather small at these high energies. Even taking into account the rather large 
error in Ci{Oi), the agreement between QCD and the data is at an impressive 
1 — 2% level. Turning now to CIPT, and repeating this procedure, but using 
a s (M T 2 ) = 0.344 ± 0.009 (A = 381 ± 16 MeV), as determined using CIPT [IT], 
we find the results shown in Fig. 4 (no vacuum condensates), and in Fig. 5 
(with condensate values as in FOPT). The experimental value of / ff is again 
reproduced within 1%, and results from both methods, FOPT and CIPT, are 
basically identical within errors. Hence, in the sequel we will only show results 
using CIPT. 

In order to study the presence of a dimension d = 2 term we shall follow 
two strategies. First, we allow for changes in the resulting pion decay con- 
stant of the order of, say, 1 MeV, i.e. a 1 % variation, and introduce a non- 
zero 02(02) with A and 04(04) fixed. This leads to the maximum value 
C 2 (0 2 ) 0.8 x 10~ 4 GeV 2 , well consistent with zero. Second, we allow 04(04} 
to change with A fixed, and introduce a dimension d = 2 term in the FESR con- 
straining the result to reproduce the experimental value of f„. In this case 
C 2 (0 2 ) remains definitely small, e.g. for C 4 (0 4 ) = 0.1 GeV 4 , and A = 397MeV 
(q s (M 2 ) = 0.353) we find C 2 (0 2 ) ^ -0.03 GeV 2 . Reducing A and C 4 (0 4 ) 
can change the sign of the dimension d = 2 term, which can reach a maxi- 
mum value C 2 (0 2 ) ^ 0.005 GeV 2 , for A = 365 MeV (a s (M 2 ) = 0.335), and 
04(04} ~ 0.05 GeV 4 . The presence of this term in the FESR turns out to be 
detrimental to the stability region, which tends to become narrower as it ap- 
proaches the experimental value of / ff at much higher energies. Our result can 
be compared with a model dependent estimate |19j . of substantially different 
order of mas nitude, C 2 (0 2 } ^ + (0.17-0.42) GeV 2 . Clearly, while our analysis 
does not categorically rule out a dimension d = 2 term in the OPE, the results 
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Figure 2: Results for f v from the FESR in the axial-vector channel, Eq. (13), 
with no condensates, and using FOPT with A = 343 MeV (a a (M?) = 0.322). 
The straight line is the experimental value of f n . 

are consistent with earlier conclusions that there is no compelling evidence for 
such a term [15], [31]- [H]. 

Next, we turn to the vector channel where the FESR is obtained from Eq.(13) 
by setting the pion pole contribution to zero. Lacking a driving term that sets 
the scale, one could analyse the FESR in several ways. We have found that 
a transparent way is to simply compare the right hand side with zero. The 
result is shown in Fig.6 for A = 397MeV (a s (M?) = 0.353), C 2 (0 2 ) = 0, 
and 6*4(04) = 0.1 GeV 4 ; the vertical axis, F(sq), is the right hand side of the 
FESR, Eq.(13). The agreement is very good in an unusually broad stability 
region. Changes in A and the corresponding dimension d = 4 condensate are 
fully compatible with the previous results from the axial-vector channel. Also, 
the conclusion about 6*2(02) being very small remains unchanged. Hence, the 
compatibility is established between results from the vector and the axial- vector 
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Figure 3: Results for f v from the FESR in the axial-vector channel, Eq. (13), 
with C 2 (0 2 ) = 0, Ca{Oa) = 0.05 GeV 4 and using FOPT with A = 343 MeV 
(a s (M%) = 0.322). The straight line is the experimental value of f„. 



channel, and in turn, their compatibility with PQCD plus no dimension d = 2 
term in the OPE, and a dimension d — 4 condensate in line with previous inde- 
pendent determinations using different weighted FESR |13|-|16]. 



We comment now on the sensitivity of the results to the choice of the cutoff 
energy s\ which enters in the integration kernel -P(s), Eq.(ll). The experimental 
errors increase appreciably as one approaches the kinematical end-point s = M 2 . 
Hence, we have chosen si high enough to include as much data as possible, but 
somewhat below the end-point so that the experimental uncertainties remain 
under control. In any case, varying s\ in the reasonable range si ~ 2.6—3.0 GeV 2 
we find that the resulting / ff changes (upwards or downwards) only by some 
0.4 — 1.3 % in the region of stability. The latter remains essentially unchanged 
in quality. 
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Figure 4: Results for f„ from the FESR in the axial- vector channel, Eq. (13), 
with no condensates and using CIPT with A = 380 Me V (a a (M?) = 0.343). 
The straight line is the experimental value of f n . 

Finally, we discuss the error A(s ) in the FESR, Eq.(13), due to the fact that 
the data might, in all probability, not be constant above the kinematical end- 
point. Assuming that PQCD sets in for s > si, we would arrive at the following 
estimate of this error 



AM _ _L {, , hM _ _ mM _ 2 M£ 1 (g_|*M j . (14) 

Numerically, we obtain A(s ) < 5 x 10~ 5 GeV 2 inside the stability region, and 
for A = 365 — 380, MeV. This is to be compared with the right hand side of 
the FESR, Eq. (13), i.e. 2 / 1 ~ 1.7 x 10~ 2 GcV 2 . This conclusion is basically 
insensitive to reasonable changes in s\, due to the fact that A (so) is some two 
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Figure 5: Results for /„■ from the FESR in the axial-vector channel, Eq. (13), 
with C 2 {0 2 ) = 0, Ci(Oi) = 0.05 GeV 2 , , and using CIPT with A = 365 MeV 
(a s (M%) = 0.335). The straight line is the experimental value of f„. 



orders of magnitude smaller than 2f%. Alternatively, relaxing the assumption 
that PQCD sets in for s > Si, we have used e + e~ data [22] in the region si — So, 
with sq in the interval s = (sj — 4) GeV 2 , in order to compute the error 
(deviation from zero) in the vector channel, i.e. 

1 f s ° 1 
A ( S °^ V = 4^J dsP(s)-ImU(s)\ e+e - DATA . (15) 

The result is A(so)|v — (10 -5 — 10 _4 )GeV 2 , i.e. one to two orders of magnitude 
smaller than F(sq) (see Fig. 6). One should notice that for sq > 4 GeV 2 the 
e + e~ data is essentially constant, albeit with very large errors, all the way up 
till the J/tp region [25] in which case A(so)|y = . 
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Figure 6: F(s ) is the right hand side of the FESR, Eq.(13), in the vector 
channel compared with zero, using CIPT with A = 397 McV (a a (M?) = 0.353), 
C 2 {0 2 ) = 0, and C 4 (e> 4 ) = 0.05 GeV 4 . 

4 Conclusions 



In this paper we have introduced a new QCD-FESR involving an integration 
kernel designed to suppress the hadronic contribution beyond the kinematical 
end-point of r-decay, i.e. the region where there is no longer data. While various 
different kinds of integration kernels in FESR have been used regularly in the 
past, this particular one is unique in its purpose as well as in the results it 
produces. In fact, this FESR allows to test the consistency between QCD and 
the r-decay vector and axial-vector hadronic spectral functions. We have found 
an excellent consistency in a remarkably wide range of energies, to wit. In the 
axial-vector channel the FESR reproduces, within errors, the experimental value 
of the pion decay constant in the broad region sq — 4 — 10 GeV 2 , for current 
values of the QCD strong coupling, and the vacuum condensates. In the vector 
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channel, where the FESR is confronted with zero, agreement is found roughly 
in the same stability region. This consistency would not be at all obvious if 
non-weighted FESR were to be employed. In fact, as it follows from Fig.l, in 
this case the saturation of the axial-vector FESR is rather poor. This is one of 
the central conclusions from this analysis. And, of course, the other is that we 
have shown that by using a suitable integration kernel it is possible to extend 
a QCD FESR analysis well beyond the kinematical end-point of data. This 
feature could become very useful in instances where the experimental data is 
only known in a narrow energy region. 
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